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Abstract 

Motivated by first-principles results for jellium and by surface-barrier shapes that are typically 
used in electron spectroscopies, the bias voltage in ballistic vacuum tunneling is treated in a heuris- 
tic manner. The presented approach leads in particular to a parameterization of the tunnel-barrier 
shape, while retaining a first-principles description of the electrodes. The proposed tunnel barriers 
are applied to Co(0001) planar tunnel junctions. Besides discussing main aspects of the present 
scheme, we focus in particular on the absence of the zero-bias anomaly in vacuum tunneling. 

PACS numbers: 72.25.Mk, 73.40. Gk, 75.47.Jn 
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I. INTRODUCTION 



At present, extensive efforts are undertaken to employ the electronic spin in 'magneto- 
electronic' devices. This aim challenges especially applied physics, but one is also con- 
cerned with model systems of spin-dependent transport in order to understand the basic 
phenomena.— Prototypical devices for studies of ballistic tunneling are planar tunnel junc- 
tions (PTJs), which consist of two magnetic electrodes separated by an insulating spacer. 
Of particular interest are the dependencies of the tunnel magneto-resistance (TMR) on the 
electronic structure of the leads and the spacer, on the width of the spacer, and on the bias 
voltage. 

The conductance of a PTJ depends on the density of states (DOS) of the electrodes and 
of the tunneling probability of the scattering channels.— The TMR can then be related to 
the spin polarization of the ferromagnetic electrodes.— Biasing, which can be viewed as a 
shift of the chemical potential of one electrode relative to that of the other, enlarges the 
range of energies in which electrons can tunnel through the spacer and introduces an energy 
dependence of the electrode spin polarization. 

State-of-the-art calculations for spin-dependent tunneling are based on the very successful 
density-functional theory (DFT).— A bias voltage, however, leads to a non-equilibrium state, 
which makes it difficult to apply DFT. An appropriate theoretical description of such a sys- 
tem would require non-equilibrium Green functions (see, e. g., Ref. |5|). Therefore, a question 
arises how one can maintain the ab-initio framework of electronic-structure calculations, in 
particular for the leads, but treat the bias voltage in a feasible manner. 

Focusing on spin-dependent ballistic tunneling through PTJs with finite bias, we investi- 
gate in the present work as a simple case tunneling through a vacuum barrier. The electronic 
properties of the electrodes from the spacer can still be computed within spin-polarized DFT. 
The crucial point is the electrostatic potential in the spacer region. Guided by first-principles 
calculations for jellium& and by theoretical models for surface barriers, we construct tunnel 
barriers that show the correct asymptotical behavior for large spacer thickness. In particular, 
one of them compares well with barrier shapes obtained ab initio for jellium. The absence of 
the zero-bias anomaly (ZBA) in vacuum tunneling of Co(0001) which was recently found by 
Ding and coworkers' 7 ^ lends itself support for an application of the proposed tunnel barrier 
(for an experimental investigation of Co PTJs with an oxide barrier, see Ref. U). We note in 
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passing that the effect of interface states on vacuum tunneling in fcc-Co(OOl) was recently 
investigated theoretically . 10 i n 

The paper is organized as follows. In Section |Hj two heuristic ways of constructing a 
tunnel barrier (jHB 21 and III B 3|) are motivated. Section III CI deals with computational as- 
pects of calculations for ballistic tunneling. Results for vacuum tunneling between Co(0001) 
electrodes are discussed in Section ITTT1 



II. THEORETICAL 



A. Surface-barrier shapes of metals 

The shape of the surface barrier of a metal was investigated in a vast amount of publi- 
cations. The possibility to calculate accurately reflected intensities in low-energy-electron 
diffraction (LEED), which is a particular surface-sensitive spectroscopy, led to several barrier 
models. Especially at very low energies (VLEED), the shape of the surface barrier has a con- 
siderable effect on the LEED /(V) spectra.— The free parameters that enter its functional 
description are fixed by fitting theoretical to experimental data, e. g., to VLEED intensities 
or to energies of surface and image-potential states.— The latter can be accessed by inverse 
or by two-photon photoelectron spectroscopy— Note that electronic-structure calculations 
using the local-density approximation (LDA) do not reproduce the correct image potential 
in the vacuum. 

Regarding electron diffraction, the classical electrostatic potential at a metal surface, with 
asymptotics V(z) ~ l/(4z), was first investigated by Mac Coll . 15 ' 16 To avoid the divergence 
at the metal surface, Cutler and Gibbons^ proposed a model potential which interpolated 
between the (constant) inner potential U of the metal and the image-charge potential in the 
vacuum region. Among the various proposed models, two became the most popular: the 
so-called JJJ barrier, named after the inventors Jones, Jennings, and Jepsen^ (see III B 31 
below), and the Rundgren-Malmstrom (RM) barrie r 19 ! 20 (for a discussion of JJJ and RM 
barriers, see Ref.Q). However, electron scattering is a dynamical process and these static 
shapes hold in principle only for the energy of interest. An energy-dependent generalization 
of the JJJ barrier suggested by Tamura and Feder proved to be successful in describing the 
image states at the Pd(llO) surface.— Further, the atomic structure at the surface leads to a 
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FIG. 1: Method of image charges for two metallic leads (grey areas). The charge q = +1 (black 
circle) is located between the left lead C (with surface at zc) and the right lead 7Z (with surface at 
z-ji). Two series of image charges are obtained by reflection at the surfaces, starting with reflection 
either at zc (blue circles) or at zn (red circles). Each image charge is indicated by the order of 
reflection (1, 2,3,.. .) and the sign of the charge (±). Only the first three orders are depicted. 

corrugated (three-dimensional) surface potential . 22 i 2 ?i 24 i 2 ^ 2 ^ However, for most applications 
the laterally (one-dimensional) and energetically invariant shapes appear to be sufficient. 

B. Construction of the tunnel barrier 

In this Section, we propose two methods of constructing the tunnel barrier. The first 
method uses the electrostatic potential of a charge between two electrodes (jHB 2|) . the 
second approach is a simple superposition of surface potentials pi B 3|) . 

1. Electrostatic potential between two metal surfaces 

Consider a planar tunnel junction with the lead L occupying the half-space ] — oo, zc], 
whereas the lead 1Z fills [zr,, oo[, with zc < z n (Fig. |TJ). The electrostatic potential V es of a 
charge q = 1 between the two semi-infinite metals can easily be obtained by the method of 
image charges.— Because each metal surface acts as mirror, one has to sum up two infinite 
series of image-charge potentials (blue and red circles in Fig. This procedure results in 

ze)zc,z n [. (1) 

Here, 7 ~ 0.577216 is Euler's constant and \I/ denotes the Digamma function.— The latter 
is the logarithmic derivative of the Gamma function T(z), ^f(z) = dlnT(z)/dz, with ^(1) = 
—7, ty(z) oc In 2; for z — > 00, and ^(z) oc —7 — 1/z for z — > + . Obviously, V cs diverges 
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for z — > zc and z — > z-r. It shows further the well-known asymptotics for the presence of a 
single metal. For example, expanding V^ s in a power series around z — z-r, yields 

1 



lim V P Jz) 



z Gj - oo, z-r . 



(2) 



ZC^-OO 4:(ZR — Z) ' 

Considering only the first-order approximation for V es , i.e., the direct images of the charge 
q (labeled 1 in Fig. [TJ), 
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(3) 



one sees that V es represents a higher barrier than Vep because the images of even order 
produce an additional repulsion ('+' in Fig. [TJ) . 



2. Image- charge potential as tunnel barrier 

The electrostatic potential between two semi-infinite jellium metals including a bias volt- 
age was calculated self-consistently within DFT by Lang.— He found that even with finite bias 
the potential in the electrodes is constant a few Bohr radii apart from the respective surfaces. 
Further, the divergence of the classical image-charge potential V^ s at zc and z-r is bridged 
over by a smooth interpolating function which shows the form of a typical LEED-motivated 
surface barrier (cf. Ill A|) . And last, application of a bias voltage apparently produces a lin- 
ear potential drop in the spacer region (cf. Fig. 2b in Ref. M). Guided by these findings we 
construct in the following a tunnel barrier by means of the classical electrostatic potential 
[eq. (JJJ)] and by LEED-type surface potentials. 

To avoid the divergences of the electrostatic potential V es , a smooth continuous interpo- 
lating function between the image-charge potential and the inner potentials of the leads, Uc 
and U-R, is used (The vacuum energy is taken as energy zero). For this paper we choose a 
Lorentzian shaped but any other reasonable shape can be used, too fsee III Aj) . The interface 
potential Vx then reads 

-U c ze]-oo,z c c ) 
ac[l + f3c(z-z c ) 2 }- 1 +lc ze[z c ,z c ] 
V,r! :) I V(z) z G [z£, z£\ ■ (4) 

a n [1 + fa(z - z^fY 1 + ln z G [z n , z n ] 
-U-r z G [z n ,oo[ 
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FIG. 2: Image-charge potential as tunnel barrier. The barrier shape is defined by the parameters 
as indicated [cf. eq. (jlj)]. The Lorentzian shapes extend over the ranges [z c , z v c \ and [z^,z^] and 
connect smoothly to the electrostatic potential in [z]* , z^\ . The latter comprises the image-charge 
potential V es [eq. ©] as well as the bias potential Vb [eq. ©]■ The inner potentials Uc = 13.61 eV 
and Un = 12.61 eV for the left and the right lead, resp., determine the bias voltage E\> to +1 eV. 

The coordinates z c c , zc , and z w c specify the positions of the onset of the interpolating 
Lorentzian, of the divergence of the potential V, and of the transition to V with respect 
to C (Fig. EJ). They have to be obtained by comparing theoretical results with other data, 
e.g., surface-state energies, VLEED spectra, etc, for the surface system (i.e., in the limit 
zqi — > oo). The parameters ace, fie, and 7£ are fixed by the conditions of smooth continuity 
in z c c and z\. Analogous considerations apply for lead 1Z. The potential V in the interior 
of the spacer can be chosen to incorporate the electrostatic potential between two metal 
electrodes, V es , and the bias voltage as well, as being discussed in the following. 

Bringing two metals so close that electrons can tunnel from metal to the other aligns the 
Fermi levels of the two leads. This energy shift is given by the contact potential $£ — $7^, 
i. e., the difference of the work functions $£ of £ and $7^ of 1Z. Note that the alignment of 
the Fermi levels is accompanied by a shift of the inner potentials, that is, e.g., Un of the 
semi-infinite system is replaced by U-n — + $7^. Since V(z) was not specified explicitly 
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FIG. 3: Dependence of the interface potential V\{ [eq. (jl))] on the lead separation. The barriers 
with electrostatic potential V es [eq. (^Q)] are shown for separations of 10, 15, 20, and 25 Bohr radii 
(blue) at zero bias (E^ = eV). In addition, a barrier with first-order approximation potential VeP 
[eq. ©] is shown for 15 Bohr radii separation (red). The inner potentials of the leads are equal 
{U c = U n = 13.61 eV). 



in eq. (jlj), it can account for the contact potential and the bias voltage, 

V(z) = V es (z) + V h (z). 



(5) 



Here, V es is the electrostatic potential from eq. (0) and H> is the bias, for which a linear drop 
over the interface region is assumed: 
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(6) 



E h z e [z^, oo[ 

This ansatz is motivated by the fact that the electric field is well screened within the elec- 
trodes but unscreened within the vacuum spacer. 



Figure H3 presents a series of tunnel barriers in dependence of the lead separation 



The heights of the barriers increase with separation (for an experimental estimation of 
the barrier height vs distance, see Ref. 3). As already mentioned, taking the first-order 
approximation instead of V cs leads to a reduced height (cf. the red line for 15 Bohr 
radii lead separation). The shape dependence on the bias is addressed in Fig. HJ For 
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FIG. 4: Dependence of the interface potential Va [eq. (JU)] on the bias. For a lead separation of 15 
Bohr radii, the right lead 1Z is biased from —6 eV to +6 eV (as indicated on the right; alternating 
blue and red lines; Uc = Un = 13.61 eV). 

rather large bias, the linear potential drop in the interface region can be clearly retrieved. 
We note in passing that the present construction produces potential shapes that compare 
qualitatively well with those obtained ab initio for jellium by Lang, 6 Further, a similar 
approach was recently used to explain the Stark shifts of surface states in scanning tunneling 
spectroscopy.^ 

One advantage of the present approach is that the height of the tunnel barrier is auto- 
matically adjusted in dependence on the lead separation and on the bias voltage. Further, 
the barrier shape shows the correct image potential asymptotics for large lead separation [cf. 
eq. In turn, the approach should not be applied for too small separations because the 
barrier shape would significantly differ from the interface potential which would be obtained 
from a self-consistent calculation for a narrow tunnel junction. This, however, could possibly 
be compensated by adjusting the parameters z^, . . . , z^ not for the semi-infinite system but 
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for the narrow junction. 



3. Superposition of surface barriers 

For large lead separations and small bias voltages, the probability of electrons to tunnel 
from one lead to the other is very small. Hence, in a self-consistent calculation for a tunnel 
junction, the tunnel barrier appears to be almost exclusively determined by the electron 
density of the respective lead and not significantly influenced by that of the other lead. This 
consideration might lead one to construct a tunnel barrier by superposition of the respective 
surface barriers, 

V i{ (z) = V c (z) + V 1z (z). (7) 

where Vc and Vr, are the surface potentials of the respective leads. Taking JJJ barriers,— 
one arrives at 

{■jt-^- — T (l — exv>\Xr(zr — z)}) z £ \zr, oof 
~a cCX p[Pctzc-z)]+l Ze}-00,Z C } 

for the surface barrier of C The values of otc and (3c are determined by requiring smooth 
continuity at z = zc- Because Uc is known from the self-consistent calculation for the surface 
system, zc and Xc remain as the only parameters to be adjusted. For the surface potential 
of 7Z one obtains an analogous form. 

The simple superposition of surface barriers appears to be problematic for hetero- 
junctions or biased junctions. In both cases, the relative energy shift of one electrode, 
say 7Z, results in a finite potential which extends into the entire other electrode C. This is 
due to the fact that the surface potential of TZ extends infinitely far into C One way to 
overcome this problem is to take the bias only as an energy shift in the interior of TZ, that 
is, to replace the inner potential U-ji by U-ji — E h . 

Figure El shows a superposition of JJJ surface barriers. The inner potential of TZ is shifted 
(cf. the arrow) by = 3 eV. Apparently, the barrier shape does not change significantly 
with bias in [zc, z-%], in contrast to the former construction (Fig.HJ). In particular, the linear 
potential drop is not observed. 
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FIG. 5: Formation of an interface barrier by superposition of surface barriers. Without bias, two 
surface barriers of JJJ type (C green and 1Z blue) are superposed to yield the interface barrier 
(black). The inner potentials are equal (Uc = Un = 13.61 eV). The application of a bias shifts the 
inner potential of 1Z (red; E\, = 3 eV, cf. the arrow) and results in the other tunnel barrier (also 
black). The positions of the image potential divergences are zc = — 7.5 ao and z-ji = 7.5 ao (Bohr 
radii), respectively (cf. the vertical dotted lines), Xc = = 1.25/ao [eq. ©]• The vacuum level 
is taken as energy zero. 

4- Resume 

The two construction recipes result in tunneling barriers with different features. While 
the more elaborate one [eq. (@J)] produces barrier shapes which are qualitatively close to 
those obtained from first-principles for jellium,— the barriers of the superposition approach 
[eq.(J7J)] lack most of these important features. In particular, the linear potential drop in the 
spacer region is missing. Both approaches can easily be extended to energy- dependent and 
corrugated (three-dimensional) tunnel barriers. 

Wang and coworkers obtained the shape of the tunnel barrier by matching two surface 
systems that were calculated for equal but opposite shifts of the work functions. 11 The 
resulting bias-dependent barriers agree well with that shown in Fig. 0J 
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C. Computational aspects of ballistic tunneling 



For the ballistic-tunneling calculations we applied the layer-KKR (Korringa-Kohn- 
Rostoker) approach of MacLaren and coworkers^ which is based on the Landauer-Biittiker 
result for the tunnel conductance.— At a given energy E and in-plane crystal momentum k\\, 
one computes the Bloch states nc and m-ji of the electrodes £ and TZ and classifies them with 
respect to their propagation direction: to the right (+) or to the left (— ). The scattering 
matrix S of the spacer S is first computed in a plane-wave basis using LEED algorithms 



(like layer-doubling and layer-stacking; see for example Ref. l31[ ) and subsequently expressed 

— * 

in terms of the scattering channels, i. e., in the Bloch-state basis. The transmission T(E t , k\\) 
at the tunnel energy E t is then a sum over all pairs of Bloch states that are incident in £ 
and outgoing in 71, 



T{E t , k\\) — fern* (Ml)*- ( 9 ) 



The tunnel conductance G(E t ) is obtained by summing over the two-dimensional Brillouin 
zone (2BZ), 

G(E t ) = G T(E t ,k\\). (10) 

fc||S2BZ 

Here, Gq = e 2 /h is the quantum of conductance which equals 2tt in atomic units. 16 Adaptive 
mesh refinement provides an efficient method to obtain accurate and well-converged 2BZ 
sums, in particular, if small parts of the 2BZ contribute significantly to the conductance.— 
With a bias voltage applied, electrons can tunnel from occupied states of one lead into 
unoccupied states of the other lead. The total conductance is then obtained by integrating 
over the energy interval given by the Fermi energies Ep of the electrodes. The averaged 
conductance thus reads 

i rmeix(E F c,E F7 z) 

G av = j- — / G(E t ) dE t . (11) 

The tunnel magneto-resistance (TMR) p is defined as the asymmetry of the (averaged) con- 
ductances for parallel (P) and antiparallel (AP) alignment of the electrode magnetizations, 

_ C av (P) - C av (AP) 
9 G av (P) + G av (AP)' 1 ] 

To treat in practice the bias voltage we proceed as follows. First, self-consistent electronic- 
structure calculations for the semi-infinite leads £ and 1Z generated muffin-tin (MT) poten- 
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tials of the bulk, of the surface, and in the vacuum region. The MT zeroes were taken 
as inner potentials Uc and Un, respectively (the MT zero is the constant potential in the 
interstitial region). For each of the leads, the MT potentials in the vacuum region were 
replaced by a smooth surface barrier. The parameters of the latter were fixed by requiring 
that the spectral densities in the surface layers were as close as possible to that of the orig- 
inal self-consistent calculation. Here, the focus laid in particular on the energy range used 
in the subsequent tunneling calculations and on surface states. For Co(0001), an important 
feature is the energy of the majority surface state at k\\ = (cf. Fig. 2b in Ref. 0; see also 
Refs.QQS). 

Having fixed the barrier parameters, the tunnel junction was built from the bulk and 
surface potentials of the two electrodes and the interface barrier [eqs. (J1J), ©> and ©]• The 
spacer S comprises all layers with potentials that differ from the respective bulk potentials. 
In example, for a Co(0001) tunnel junction the first four layers on either side of the smooth 
tunnel barrier were used. The bias was taken into account by shifting the inner potential of 
one of the leads (muffin-tin zero) and determining the barrier shape [eq. (JTj)]. The smooth 
barrier Va was treated as a single layer in the multiple-scattering calculations. Its scattering 
matrix S was obtained within the propagator formalism.— 

As usual for the KKR method, a small imaginary part r\ has to be added to the energy E, 37 
leading in general to complex wavenumbers fcj_ . 38 i 39 Therefore, the electrode eigenfunctions 
are no longer true Bloch states but become evanescent states [^(fcjj ^ 0]. Eigenstates 
stemming from Bloch states [3(A;^) = for 77 = 0] show typically the smallest Q(k±) and 
can therefore be separated from evanescent states [Q^&jJ 7^ for rj = 0].— In the spacer 
iS, the nonzero rj leads to damping in addition to the intrinsic one, artificially enhancing 
the decay of the conductance with spacer thickness. Further, the scattering matrix 5* is no 
longer unitary and, hence, the total current is not conserved. Therefore, one has to choose rj 
carefully in order to produce reliable results. We found that a value of rj = 10~ 4 eV produces 
no considerable artefacts. 

The Landauer-Biittiker approach used here avoids the computation of the Green function 
of the complete system, which is in particular problematic for a non-equilibrium system. 
Considering the asymptotic transmission channels (Bloch states), states that are localized 
at the barrier do not contribute to the transmission. 

Recently, Davis and MacLaren reported on model calculations for spin-dependent tunnel- 
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ing at finite bias.— In their work, however, the electronic structure of the Fe electrodes was 
approximated by plane waves, whereas the barrier was assumed as step-like with a linear 
drop. Although conceptual similar, our approach goes beyond that work. First, the elec- 
trodes are treated on first-principles level. Second, the barrier shows the correct asymptotics 
(for the free surfaces) and, once the shape parameters being fixed, depends automatically 
on both lead separation and bias. 



III. RESULTS FOR Co(0001) 

Recently, Ding and coworkers investigated the bias-voltage dependence of the TMR with 
a spin-polarized scanning tunneling microscope (STM).— In contrast to tunneling through 
oxide barriers, they observed no zero-bias anomaly (ZBA), i. e., a (rather) sharp maximum of 



the TMR at zero bias (see, for example, Ref. l4ll). With a vacuum barrier replacing an oxide 



barrier, the TMR appeared to be almost constant. This finding suggests that the ZBA is 
mainly due to imperfections in oxide barriers, rather than to scattering at magnons and spin 
excitations (in the leads). Further, the so-called DOS effect, i.e., the energy dependence of 
the spin-resolved density of states of the leads, proved to be small in the case of Co(0001). 

The experimental findings of Ding et al. were corroborated by ballistic tunneling calcu- 
lations for planar Co(0001) junctions as sketched in Section III CI The tunnel barrier was 
taken as a superposition of surface barriers pi B 3J) . In the present work, we focus on the 
more elaborate image-charge potential pi B 2|) . 

The transmission T(Ep, kn), eq. Q, depends on the relative orientation of the lead mag- 
netizations (P and AP), as is shown for eV bias (tunneling at Ep) in Fig. |H1 For the 
chosen lead separation of 7.52 A, only those Bloch states with a k\\ in the central part of the 
2BZ contribute significantly to the transmission. The normal component of the wavevec- 



tor within the tunnel barrier, kj_(z) = y2[E-p — Vif(z)] — fcjj , is imaginary and gives rise to 
strongly evanescent states in the tunnel barrier for Bloch states with large fen, and, thus, 
to a small transmission. For Bloch states with kn near T (i.e., k\\ = 0), the decay within 
the barrier is less and the transmission can be larger. In total, this results in a 'focusing' of 
T(E t , k\\) at the 2BZ center. 

Both the P and the AP case show minor transmission close to T. These minima are 
surrounded by ring-like structures of increased transmission. The maximum P transmission 
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FIG. 6: Transmission T(Ep, fen) of vacuum tunneling in Co(0001) for P (a, top) and AP (b, bottom) 
alignment of the lead magnetizations. For eV bias, the lead separation was chosen as 7.52 A. 
The maximum transmission is about 0.01 (P alignment. Both panels share the same scale). Note 
that only the central part of the two-dimensional Brillouin zone is displayed. 

is larger than for AP alignment (by a factor of about 10). But the AP transmission displays 
a broader ring compared to the P transmission. 

When integrated over the 2BZ, one finds that G(P) > G(AP) (cf. the black symbols in 
Fig- HI for zero bias). With increasing bias, the conductances for tunneling at E-p decrease. 
Since ballistic tunneling is a phase coherent process, shifting of the electronic states of one 
electrodes relative to those of the other by the bias, might reduce the phase coherence. 
Or the spin-dependent DOS in the relevant region of the 2BZ decreases with energy. We 
checked the spectral density carefully but found no significant feature that would corroborate 
unequivocally the latter explanation. 

The conductances increase with tunnel energy E t (cf. the data for 0.2 to 1.0 eV bias). 
This fact might be explained by the electronic structure of the leads or by a reduction of 
the effective barrier width and height (cf. Figs. Eland EJ). Indeed, the radius of the ring-like 
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FIG. 7: Ballistic conductance G (in units of Go, the quantum of conductance) vs. tunnel energy 
E t in Co(0001) for P (a, bottom) and AP (b, top) alignment. The bias E^ ranges from 0.0 eV 
(black) to 1.0 eV (orange) and is indicated at the top of each dataset. 

structure in the 2BZ which contributes most to the conductance (cf. Fig. [BJ) gets larger with 
tunnel energy, hence increasing the contributing area. This is consistent with the focusing 
effect mentioned earlier. The AP conductances in particular can be represented reasonably 
well by parabolae. 

A further interesting feature is the increase for biases of 0.6, 0.8, and 1.0 eV that occurs 
for P alignment at tunnel energies E t around 0.0, 0.2, and 0.4 eV, respectively (Fig. EK)- 
Inspection of the transmissions and of the spectral density at Ep — 0.6 eV produced no 
significant feature that would explain this behavior (This is corroborated by findings of 
LeClair et at, Ref. M). The feature occurs also for AP alignment (Fig. EJd) but not as 
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FIG. 8: Magneto-resistance of vacuum tunneling in Co(0001). (a) Averaged conductance G av (in 
units of Go, the quantum of conductance) vs. bias E\> for P (blue) and AP (red) alignment [cf. 
eq. (b) Tunnel magneto-resistance p vs. E^ [eq. (|12|)]. 

pronounced as for the P case. 

The increase with bias compensates the decrease for E t = Ep, as is shown for the averaged 
conductances G av in Fig. |BK Whereas G av (P) decreases slightly (with a small minimum at 
zero bias), G av (AP) increases with \E h \. Therefore, the resulting TMR p [eq. drops 
with bias, too. However, the decrease which is about 15 % at 0.6 eV is much less than that 
observed for oxide barriers. In the latter case, the TMR drops by 50 % to 80 % at 0.6 eV.— 
Being due to the details of the electronic structure in the Co leads, one could term the drop 
in Fig. |Hb as 'DOS effect', rather than as zero-bias anomaly. 

Since inelastic processes (scattering at magnons, spin excitations) are not included in 
our theory, one can conclude that the ZBA found in tunnel junctions with oxide barriers 
can be attributed to defect scattering in the oxide barrier. This finding is consistent with 
the fact that the ZBA decreases with the improvement of the preparation techniques for 
ferromagnet-oxide interfaces (see Ref . and references therein) . 

In a previous investigation,"^ we used the superposition approach pi B 3)1 for the tunnel 
barrier. There, both the averaged conductances and the TMR were almost constant for 
biases up to 0.5 eV. Comparing with the present results that were obtained within the 
image-potential approach (jll B 2)1 . one has to keep in mind that details of the calculations 



16 



differ (e.g., the k\\ mesh). However, these have only minor influence. The most striking 
difference is the shape of the tunnel barrier which is varied in two aspects. First, the JJJ 
barrier used in Ref. la is rather smooth with respect to the interpolating Lorentzian chosen 
here. Generally speaking, the latter produces a larger reflection. Second, the shape in the 
central part of the barrier differs. In particular for large lead separations, the barrier height 
becomes important (cf. Figs. El and EJ). Further, the linear bias potential which is missing in 
the superposition approach is expected to have a non-negligible effect. Therefore, details of 
the barrier are expected to have significant influence on the tunnel magneto-resistance. 

IV. CONCLUDING REMARKS 

Tunnel calculations provide a rather indirect test of the proposed barrier shapes. A more 
direct one would be to compare theoretical energy positions and linewidths of so-called field- 
emission resonances^ 2 - with experimental ones. These electronic states can be viewed as 
surface states that are trapped between the bulk (in the presence of a bulk-band gap) and 
the tunnel barrier between sample and an STM tip. The field-emission resonances show up 
as sharp maxima in the differential conductance and depend — like the shape of the tunnel 
barrier — on both bias voltage and tip-sample separation. 

As a possible extension of the present work, one could think of a treatment of tunnel 
junctions with 'filled' spacers (instead of vacuum), in particular with oxide barriers. Further, 
work is in progress to describe the tunneling with bias voltage fully on an ab-initio level. 
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